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Abstract. Consider a plane monochromatic wave incident on a semi-infinite 
periodic structure. What happens if the normal component of the transmitted 
wave group velocity vanishes? At first sight, zero normal component of the 
transmitted wave group velocity simply implies total reflection of the incident 
wave. But we demonstrate that total reflection is not the only possible outcome. 
Instead, the transmitted wave can appear in the form of a frozen mode with very 
large diverging amplitude and either zero, or purely tangential energy flux. The 
field amplitude in the transmitted wave can exceed that of the incident wave by 
several orders of magnitude. There are two qualitatively different kinds of frozen 
mode regime. The first one is associated with a stationary inflection point of 
electromagnetic dispersion relation. This phenomenon has been analyzed in our 
previous publications. Now, our focus is on the frozen mode regime related to a 
degenerate photonic band edge. An advantage of this new phenomenon is that it 
can occur in much simpler periodic structures. This spectacular effect is extremely 
sensitive to the frequency and direction of propagation of the incident plane wave. 
These features can be very attractive in a variety practical applications, such 
as higher harmonic generation and wave mixing, light amplification and lasing, 
highly efficient superprizms, etc. 
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1. Introduction 

Wave propagation in spatially periodic media, such as photonic crystals, can be 
qualitatively different from any uniform substance. The differences are particularly 
pronounced when the wavelength is comparable to the primitive translation L of the 
periodic structure mGHIIlEI03|I3|7]. The effects of strong spatial dispersion culminate 
when the group velocity u = duo/dk of a traveling Bloch wave vanishes. One reason for 
this is that vanishing group velocity always implies a dramatic increase in density of 
modes at the respective frequency. In addition, vanishing group velocity also implies 
certain qualitative changes in the eigenmode structure, which can be accompanied by 
some spectacular effects in electromagnetic wave propagation. A particular example 
of the kind is the frozen mode regime associated with a dramatic enhancement of the 
wave transmitted to the periodic medium [HI El QUI 021 QUI ■ There are at least 
two qualitatively different modifications of the frozen mode regime, each related to 
a specific singularity of the electromagnetic dispersion relation. Both effects can be 
explained using the simple example of a plane electromagnetic wave normally incident 
on a lossless semi-infinite periodic structure, as shown in Fig. ^ 
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Figure 1. Plane wave normally incident on a semi-infinite photonic crystal. The 
subscripts /, R, and T refer to the incident, reflected and transmitted waves, 
respectively. In all cases, the amplitude of the incident wave is unity. 

The frozen mode regime of the first kind is associated with a stationary inflection 
point on the k — u> diagram shown in Fig. Ejb) . In the vicinity of stationary inflection 
point, the relation between the frequency u> and the Bloch wave number k can be 
approximated as 

uj — cjo oc (k — fcfj) 3 ■ (1) 

A monochromatic plane wave of frequency close to ujq incident on semi-infinite 
photonic crystal is converted into the frozen mode with infinitesimal group velocity 
and dramatically enhanced amplitude, as illustrated in Fig. The saturation value 
of the frozen mode amplitude diverges as the frequency approaches its critical value 
ujq. Remarkably, the photonic crystal reflectivity at u> = ujq can be very low, implying 
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Figure 2. Schematic examples of dispersion relations displaying different 
stationary points: (a) a regular band edge (RBE), (b) a stationary inflection 
appoint (SIP), (c) a degenerate band edge (DBE). 



that the incident radiation is almost totally converted into the frozen mode with zero 
group velocity, diverging amplitude, and finite energy flux close to that of the incident 
wave [HI HOI [TT]. 

A qualitatively different kind of frozen mode regime is expected in the vicinity 
of a degenerate photonic band edge shown in Fig. EJc). This case is the main focus 
of our investigation. At frequencies just below u>d, the dispersion relation can be 
approximated as 

u>d — u> oc (fc — fc d ) 4 , at u> ^ uid- (2) 

Contrary to the case of stationary inflection point , in the vicinity of a degenerate 
band edge the photonic crystal becomes totally reflective. But at the same time, 
the steady-state field inside the periodic medium (at z > 0) develops a very large 
amplitude, diverging as the frequency approaches its critical value Wd- Such a behavior 
is illustrated in Fig. 0] The frozen mode profile below and above the degenerate band 
edge frequency ojd is different, ft has a large saturation value at frequencies located 
inside the transmission band (at uj < uJd), as seen in Fig. 0fa) and (b). This saturation 
value diverges as uj — ► ojd — 0. By contrast, at frequencies inside the band gap (at 
uj > uJd), the field amplitude initially increases dramatically with the distance z from 
the surface, but then vanishes as the distance z further increases, as seen in Fig. 0Jd 

-f). 

Figs. |3]and0]describe the frozen mode profile in hypothetical lossless semi-infinite 
periodic media. In the case of a photonic crystal with finite thickness, the frozen 
mode profile remains unchanged in the leftmost portion of the periodic structure in 
Fig. ^ In the opposite, rightmost part of the photonic crystal, the frozen mode 
amplitude vanishes, as illustrated in Fig. Additional factors limiting the frozen 
mode amplitude include structural imperfections of the periodic array, absorption, 
nonlinearity, deviation of the incident radiation from plane monochromatic wave, etc. 
Still, with all these limitations in place, the frozen mode regime can be very strong. 

Not every periodic structure can support the frozen mode regime at normal 
incidence. Generally, the physical conditions for the frozen mode regime are the 
same as the conditions for the existence of the respective stationary point or 
of the dispersion relation. In cither case, a unit cell of the periodic layered structure 
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Figure 3. Smoothed profile of the frozen mode at six different frequencies 
in the vicinity of stationary inflection point: (a) ui = ujq — 10~ 4 c/L, (b) 
uj = u> - 10~ 5 c/L, (c) u) = oj , (d) uj = u) + 10- 5 c/L, (c) u> = u> + 10" 4 c/L, 
(f) uj = ujq + 10 _3 c/L. In all cases, the incident wave has the same polarization 
and unity amplitude. The distance z from the surface of semi-infinite photonic 
crystal is expressed in units of L. Physical parameters of the periodic sturcture 
are specified in 1781 and 1791 , 



must contain at least three layers, of which two must display a misaligned in-plane 
anisotropy, as shown in Fig. The difference, though, is that a stationary inflection 
point also required the presence of magnetic layers with strong nonreciprocal 
circular birefringence |S] • No magnetic layers are needed for a degenerate band 
edge (0), which constitutes a major practical advantage of the respective frozen mode 
regime. In photonic crystals with three dimensional periodicity, the presence of 
anisotropic constitutive component may not be necessary. A detailed comparative 
analysis of the above two modifications of the frozen mode regime at normal incidence 
is carried out in the next section. The emphasis is on the physical conditions under 
which these phenomena can occur. 
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Figure 4. Smoothed profile of the frozen mode at six different frequencies in the 
vicinity of degenerate band edge: (a) ui = — 10~ 4 c/L, (b) uj = u) d — W~ e c/L, 
(c) u = ui d , (d) u> = ui d + 10~ 6 c/L, (e) ui = uj d + 10~ 5 c/L, (f ) uj = ui d + 10" 4 c/L. 
In the transmission band (at u> < uj d ), the asymptotic field value diverges as 
uj —> uj d . By contrast, in the band gap (at uj > ui d ), the asymptotic field value is 
zero. The amplitude of the incident wave at z < is unity. Physical parameters 
of the periodic structure used in computations are specified in Section 5. 

In Section 3, we turn to the case of oblique wave propagation. The frozen mode 
regime at oblique incidence can occur when the normal component of the transmitted 
wave group velocity vanishes, while its tangential component remains finite. In such a 
case, the transmitted wave is an abnormal grazing mode with a dramatically enhanced 
amplitude and nearly tangential energy flux. A significant advantage of the oblique 
modification of the frozen mode regime is that it can occur in much simpler periodic 
structures, compared to those supporting the frozen mode regime at normal incidence. 
Examples are shown in Figs. and |H1 The presence of anisotropic layers is still 
required. 

Yet another interesting modification of the frozen mode regime are abnormal 



Frozen light in photonic crystals with degenerate band edge 





Figure 5. Smoothed profile of the frozen mode in periodic layered structures 
composed of different number N of unit cells L. The frequency is equal to that of 
the degenerate band edge. The initial rate of growth of the frozen mode amplitude 
is virtually independent of N and described by 1341 . The limiting case (d) of the 
semi-infinite structure is also shown in Fig. 0Jc). In all cases, the incident wave 
has the same polarization and unity amplitude. 



subsurface wave. Such waves can exist at band gap frequencies close to a degenerate 
photonic band edge. Regular surface waves usually decay exponentially with the 
distance from the surface in either direction. By contrast, abnormal subsurface waves 
are extremely asymmetric. They do decay rapidly outside the photonic crystal. 
But inside the periodic medium, their amplitude sharply increases, and reaches its 
maximum at a certain distance from the surface. Only after that the field amplitude 
begins a slow decay, as the distance from the surface further increases. The profile 
of a subsurface wave is similar to that of the frozen mode above the degenerate band 
edge in Fig. 0Jd - f). This phenomenon is briefly discussed in section 4. 

In Section 5 we discuss the physical requirements to the spatially periodic arrays 
capable of supporting the frozen mode regime, both at normal and oblique incidence. 
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We also present a detailed description of the periodic layered structures used in our 
numerical simulations. 

Finally, in Section 6 we summarize the results and discuss some physical 
limitations of the frozen mode regime. 
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Figure 6. Periodic stack capable of supporting k — ui diagram with a DBE. A unit 
cell L includes three layers: two birefringent layers A\ and A2 with misaligned 
in-plane anisotropy, and one isotropic B layer. In order to support a DBE, the 
misalignment angle <f> between adjacent anisotropic layers Ai and A2 must be 
different from and it/2. A detailed description of this periodic structure is given 
in the Section 5. 
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Figure 7. Periodic layered structure with two layers A and B in a primitive cell 
L. The A layers (hatched) are anisotropic with one of the principle axes of the 
dielectric permittivity tensor making an oblique angle with the normal z to the 
layers (e xz ^ 0). The B layers are isotropic. The x — z plane coincides with the 
mirror plane of the stack. This structure can support axial dispersion relation 
u> (fc z ) with stationary inflection point 1441 . provided that k X; k y 7^ 0. 
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Figure 8. Periodic layered structure with two layers A and B in a unit cell L. 
The A layer has inplane anisotropy 1761 . while the B layer can be isotropic. This 
stack can display axial dispersion relation ui (fc z ) with a degenerate band edge 
1451 . provided that k x ,k y ^ 0. 



2. The physical nature of the frozen mode regime 

The essence of the frozen mode regime can be understood from the simple example 
of a plane monochromatic wave normally incident on a semi-infinite periodic layered 
structure, as shown in Fig. ^ An important requirement, though, is that some of 
the layers display a misaligned in-plane anisotropy as shown in the example in Fig. HJ1 
Below we present a comparative analysis of two different kinds of frozen mode regime. 
Although throughout this section we only consider the case of normal incidence, in 
the next section we will show that most of the results and expressions remain virtually 
unchanged in a more general case of the frozen mode regime at oblique propagation. 
One difference, though, is that at oblique incidence, the frozen mode regime can occur 
in much simpler structures. This can have a big advantage in practical terms. 

To start with, let us introduce some basic notations and definitions. Let ^r, 
and be the incident, reflected and transmitted waves, respectively. Assume for now 
that all three monochromatic waves propagate along the z axis normal to the surface 
of semi-infinite periodic layered structure in Fig. ^ Electromagnetic field both inside 
(at z > 0) and outside (at z < 0) the periodic stack is independent of the x and y 
coordinates. The transverse field components can be represented as a column- vector 



*(*) = 



E x {z) 

Ey(z) 

H x (z) 

Hy(z) 



(3) 



where E (z) and H (z) are time-harmonic electric and magnetic fields. All four 
transverse field components in © are continuous functions of z, which leads to the 
following standard boundary condition at z = 



*t (0) = */ (0) + (0) 



(4) 
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Assume also that anisotropic layers of the periodic array have an in-plane 
anisotropy (|69|) . in which case the fields E (z) and H (z) are normal to the direction 
of propagation 

E(z) J- z,H (z) Lz, (5) 

and the column vector JSJl includes all nonzero field components. Note that the 
polarizations of the incident, reflected and transmitted waves can be different, because 
some of the layers of the periodic array display an in-plane anisotropy, as shown in the 
example in Fig. [U The presence of anisotropic layers is essential for the possibility of 
frozen mode regime. 

In periodic layered media, the electromagnetic eigenmodes {z) are usually 
chosen in the Bloch form 

* fc (z + L) = e lkL V k (z) , (6) 

where the Bloch wavenumber k is defined up to a multiple of 2n/L. The 
correspondence between lo and k is referred to as the Bloch dispersion relation. Real 
k correspond to propagating (traveling) Bloch modes. Propagating modes belong to 
different spectral branches lo (k) separated by frequency gaps. The speed of a traveling 
wave in a periodic medium is determined by the group velocity 4; 

u = dio/dk. (7) 

Normally, each spectral branch lo (k) develops stationary points lo s — lo (k s ) where the 
group velocity Q of the corresponding propagating mode vanishes 

dio/dk — 0, at lo — lo s = lo (k s ) . (8) 

Examples of different stationary points are shown in Fig. [3 where each of the 
frequencies lo 9 , loq and lo^ is associated with zero group velocity of the respective 
traveling wave. Stationary points play essential role in the formation of frozen 
mode regime. 

By contrast, evanescent Bloch modes are characterized by complex wavenumbers 
k = k' + ik" . Evanescent modes decay exponentially with the distance z from the 
boundary of semi-infinite periodic structure. Therefore, under normal circumstances, 
evanescent contribution to the transmitted wave (z) can be significant only in close 
proximity of the surface. The situation can change dramatically when the frequency lo 
approaches one of the stationary point values lo s . At first sight, stationary points (JSJ 
relate only to propagating Bloch modes. But in fact, in the vicinity of every stationary 
point frequency to s , the imaginary part k" of the Bloch wavenumber of at least one of 
the evanescent modes also vanishes. As a consequence, the respective evanescent mode 
decays very slowly, and its role may extend far beyond the photonic crystal boundary. 
In addition, in the special cases of interest, the electromagnetic field distribution "J (z) 
in the coexisting evanescent and propagating eigenmodes becomes very similar, as lo 
approaches io s . This can result in spectacular resonance effects, such as the frozen 
mode regime. What exactly happens in the vicinity of a particular stationary point 
JHJ essentially depends on its character and appears to be very different in each of the 
three cases presented in Fig. |3] 

In the next subsection we present a simple qualitative picture of the frozen 
mode regime based solely on energy conservation consideration. This will allow us 
to highlight the difference between the cases of stationary inflection point and 
degenerate band edge @. Then, we discuss the physical nature of the frozen mode 
regime. 
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2.1. Energy density and energy flux at frozen mode regime 

Let Si, Sr and St be the energy flux in the incident, reflected and transmitted waves 
in Fig. n The transmission and reflection coefficients of a lossless semi-infinite medium 
are defined as 

St S r 

where 

Si + Sr = St, p = X — t. 

With certain reservations, the energy flux St of the transmitted travelling wave can 
be expressed as 

S T = W T u, (10) 

where u is the group velocity, which is also the energy velocity, and Wt is the energy 
density 

W T oc |*t| 2 ■ 

Evanescent modes do not contribute to the normal energy flux St in the case of 
a lossless semi-infinite periodic structure. Besides, evanescent contribution to the 
transmitted wave becomes negligible at a certain distance z from the surface. The 
simple expression (|1UH may not apply when the transmitted wave involves two or 
more propagating Bloch modes, but we will not deal with such a situation here. 

Vanishing group velocity u in (|10|l implies that the transmitted wave energy flux 
St also vanishes, along with the respective transmission coefficient t in ©. The 
only exception could be if the energy density Wt of the transmitted wave increases 
dramatically in the vicinity of the stationary point frequency. In other words, if Wt 
in IjlOJI grows fast enough, as to approaches co s , the product Wtu in H10[) can remain 
finite even at to = lu s . In such a case, a significant fraction of the incident radiation 
can be converted into the slow mode inside the semi-infinite periodic array. The effect 
of a dramatic growth of the transmitted wave amplitude in the vicinity of a stationary 
point (JSJ will be referred to as the frozen mode regime. The possibility of such an 
effect is directly related to the character of a particular stationary point. From this 
point of view, let us consider three different situation presented in Fig. [3 



2.1.1. Regular band edge We start with the simplest case of a regular photonic band 
edge (RBE) in Fig. |2Ia) . It can be found in any periodic array, including any periodic 
layered structure. Just below the band edge frequency w 9 , the dispersion relation can 
be approximated by a quadratic parabola 

uj g — co cx (k — kg) 2 , at w £ to g . (11) 

This yields the following frequency dependence of the propagating mode group velocity 
u inside the transmission band 

u = — oc (kg — k) cx (uj g — lu) 1 ^ 2 , at uj iS uj g . (12) 

(IK 

Due to the boundary condition |0J, the amplitude of the transmitted propagating 
Bloch mode remains finite and comparable to that of the incident wave. Therefore, the 



Frozen light in photonic crystals with degenerate band edge 



11 



energy flux (|10|) associated with the transmitted slow mode vanishes, as uj approaches 

UJg 

c w f (wo-w) 1 ' 2 , at uj uj„ , 1q ^ 
S T = W T u « <, 9 „ ' ' ~ 9 . (13) 
l_ 0, at uj > ujg 

Formula (| 1 31) expresses the well-known fact that in the vicinity of a regular photonic 
band edge, a lossless semi-infinite photonic crystal becomes totally reflective. 

2.1.2. Stationary inflection point A completely different situation occurs in the 
vicinity of a stationary inflection point in Fig. |2fb). At normal propagation, such 
a point can be found in periodic layered structures involving anisotropic and magnetic 
layers jSJ|U], as well as in some photonic crystals with 2- and 3-dimensional periodicity. 
In the vicinity of a stationary inflection point ujq, the dispersion relation can be 
approximated by a cubic parabola The propagating mode group velocity u 
vanishes as uj approaches ujo from either direction 

u = — oc (k — ko) 2 oc {uj — ciJo) 2 ^ 3 ■ (14) 
dk 

But remarkably, the amplitude of the transmitted propagating mode increases so that 
the respective energy density Wt diverges as uj — > ujq 

Wt oc (uj - uj y 2/3 ■ (15) 

The expression l|14ll together with (|15|l yield that the energy flux of the transmitted 
slow mode remains finite even at uj = ujq 

St = Wtu ~ Si, at uj w ujq. (16) 

The latter implies that the incident light is converted to the frozen mode with 
infinitesimal group velocity (|14|) and diverging amplitude l|15(l . This result was first 
reported in UJ • A consistent analytical description of the asymptotic behavior of the 
transmitted field amplitude in the vicinity of a stationary inflection point was carried 
out in [T3|. 

2.1.3. Degenerate band edge Let us turn to the case of a degenerate band edge in 
Fig. Efc). At normal propagation, such a point can be found in dispersion relation of 
periodic layered structures with misaligned anisotropic layers. An example is shown 
in Fig. |B| Just below the degenerate band edge frequency ujd, the dispersion relation 
uj (k) can be approximated by a biquadratic parabola (J2J • This yields the following 
frequency dependence of the propagating mode group velocity inside the transmission 
band 

u = — tx (k d - kf cx (uj d - uj) 3/4 , at uj < uj d . (17) 
dk 

Analysis shows that the amplitude of the transmitted slow mode in this case diverges, 
as the frequency approaches the band edge value 

W T oc|w d - w\ , atw£w flBE , (18) 

which constitutes the frozen mode regime. But the energy density l|18|l does not grow 
fast enough to offset the vanishing group velocity i(TT|l . The expressions (fTTjl and (|TH|l 
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together with (fTUfl yield for the energy flux 

St = W T u oc I y a > ' ~ , 19) 
L 0, at uj > uj d 

implying that, in spite of the diverging energy density JTSJ, the energy flux of the 
transmitted slow wave vanishes, as uj approaches Ud- 

The situation at a degenerate band edge J2J can be viewed as intermediate 
between the frozen mode regime at a stationary inflection point |JT]|. and the vicinity 
of a regular band edge i|llfl . Indeed, on the one hand, the incident wave at u> = uj d 
is totally reflected back to space, as is the case at a regular band edge. On the other 
hand, the transmitted field amplitude inside the periodic medium diverges as uj -> w^, 
which is similar to what occurs at a stationary inflection point. 

The above consideration does not explain the nature of the frozen mode regime, 
nor does it address the problem of the Bloch composition of the frozen mode. These 
questions are the subject of the next subsection. 

2.2. Bloch composition of frozen mode 

In a periodic layered structure, at any given frequency uj, there are four electromagnetic 
eigenmodes with different polarizations and wavenumbers. But in the setting of 
Fig. ^ where the semi-infinite periodic array occupies the half-space z > 0, the 
transmitted wave is a superposition of only two of the four Bloch eigenmodes. Indeed, 
neither the propagating modes with negative group velocity, nor evanescent modes 
exponentially growing with the distance z from the surface, contribute to {z) in 
this case. Generally, one can distinguish three different possibilities. 

(i) Both Bloch components of the transmitted wave are propagating modes 

* T (*) = * prl (*) + % r2 (z) , z>0. (20) 

^pri (z) and ^ P r2 (z) are two propagating Bloch modes with different real 
wavenumbers k\ and k 2 and different group velocities iti > and u 2 > 0. This 
constitutes the phenomenon of double refraction, provided that U\ and u 2 are 
different. The other two Bloch components of the same frequency have negative 
group velocities and cannot contribute to the transmitted wave ^t- 

(ii) Both Bloch components of <3/t are evanescent 

* T (*) = (z) + *ev2 (Z) , Z>0- (21) 

The respective two values of k are complex with positive imaginary parts k" > 0. 
This is the case when the frequency uj falls into photonic band gap at uj > uj g 
in Fig. Ufa) or at uj > ujd in Fig. |2Jc). The fact that k" > implies that the 
wave amplitude decays with the distance z from the surface. In the case (|21|) . 
the incident wave is totally reflected back to space by the semi-infinite periodic 
structure. 

(iii) One of the Bloch components of the transmitted wave is a propagating mode 
with u > 0, while the other is an evanescent mode with k" > 



* T (z) = ly (z) + y ev (z), z> 0. 



(22) 
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For example, this is the case at ui ~ ujq in Fig. Gib), as well as at to < tv g in Fig. 
ffia) and at u> < uj^ in Fig. Hfc). As the distance 2 from the surface increases, 
the evanescent contribution $ CT in IL'l't decays as exp(— zk"), and the resulting 
transmitted wave (-2) turns into a single propagating Bloch mode ^ pr - 

Propagating modes with u > and evanescent modes with k" > are referred 
to as forward waves. Only forward modes contribute to the transmitted wave (z) 
in the case of a periodic semi-infinite stack. The propagating modes with u < and 
evanescent modes with k" < are referred to as backward waves. The backward waves 
never contribute to the transmitted wave inside the periodic semi-infinite stack in 
Fig. n This statement is based on the following two assumptions: 

- The transmitted wave 'J't and the reflected wave are originated from the 
plane wave 'J/ incident on the semi- infinite photonic slab from the left, as shown 
in Fig. [U 

- The layered array in Fig. ^occupies the entire half-space and is perfectly periodic 
at z > 0. 

If either of the above conditions is violated, the electromagnetic field inside the 
periodic stack can be a superposition of four Bloch eigenmodes with either sign of 
the group velocity u of propagating contributions, or either sign of k" of evanescent 
contributions. This would be the case if the periodic layered array in Fig. ^had some 
kind of structural defects or a finite thickness. At the end of this section we briefly 
discuss how it would affect the frozen mode regime. 

Note also that the assumption that the transmitted wave (z) is a superposition 
of propagating and/or evanescent Bloch eigenmodes may not apply if the frequency 
u exactly coincides with one of the stationary point frequencies JSJ. For example, at 
frequency uiq of stationary inflection point , there are no evanescent solutions to the 
Maxwell equations {52} , and the transmitted wave (z) is a (non-Bloch) Floquet 
eigenmode linearly growing with z Elj ■ Similar situation occurs at frequency ujd 
of degenerate band edge J2J). The term "non-Bloch" means that the respective field 
distribution does not comply with the relation ©. At the same time, at any general 
frequency, including the vicinity of any stationary point ©, the transmitted wave 
i&T (z) is a superposition of two Bloch eigenmodes, each of which is either propagating, 
or evanescent. 

In all three cases (|2()l-l2^|l. the contribution of a particular Bloch eigenmode to 
the transmitted wave depends on the polarization "J/ of the incident wave. One 
can always choose some special incident wave polarization, such that only one of the 
two forward Bloch modes is excited and the transmitted wave is a single Bloch 
eigenmode. In the next subsection we will see that there is no frozen mode regime in 
the case of a single mode excitation. This fact relates to the very nature of the frozen 
mode regime. 

Knowing the Bloch composition of the transmitted wave we can give a semi- 
qualitative description of what happens when the frequency u> of the incident wave 
approaches one of the stationary points {BJ in Fig. More consistent analysis based 
on the Maxwell equations is outlined in Section 4. 

2.2.1. Regular photonic band edge We start with the simplest case of a regular 
photonic band edge. There are two different possibilities in this case, but none of 
them is associated with the frozen mode regime. The first one relates to the trivial 
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case where none of the layers of the periodic structure displays an in-plane anisotropy 
or gyrotropy. As the result, all eigenmodes are doubly degenerate with respect to 
polarization. A detailed description of this case can be found in the extensive literature 
on optics of stratified media El ■ Slightly different scenario occurs if some of the 
layers are anisotropic or gyrotropic and, as a result, the polarization degeneracy is 
lifted. Just below the band edge frequency ui g in Fig. EJa), the transmitted field 
(z) is a superposition l)22|l of one propagating and one evanescent Bloch modes. 
Due to the boundary condition the amplitude of the transmitted wave at z = is 
comparable to that of the incident wave. In the case of a generic polarization of the 
incident light, the amplitudes of the propagating and evanescent Bloch components 
at z = are also comparable to each other and to the amplitude of the incident light 

|*pr (0)| ~ |*„, (0)| ~ |*/| , at to < UJ g . (23) 

As the distance z from the surface increases, the evanescent component ^ ev (z) 
decays rapidly, while the amplitude of the propagating component remains constant. 
Eventually, at a certain distance from the slab surface, the transmitted wave (z) 
becomes very close to the propagating mode 

* T (z) w ^ pr (z), at z » L, oj < uj g . (24) 

The evanescent component \I/ ev of the transmitted wave does not display any 
singularity at the band edge frequency to g . The propagating mode ^5> pr does develop a 
singularity associated with vanishing group velocity at uj — > uj g — 0, but its amplitude 
remains finite and comparable to that of the incident wave. At w > u) g , this 
propagating mode turns into another evanescent mode in l|21|) . The bottom line is that 
none of the Bloch components of the transmitted wave develops a large amplitude in 
the vicinity of a regular photonic band edge. There is no frozen mode regime in this 
trivial case. 

2.2.2. Stationary inflection point A completely different situation develops in the 
vicinity of a stationary inflection point Q of the dispersion relation. At lo « ujq, the 
transmitted wave is a superposition (1221) of one propagating and one evanescent 
Bloch component. In contrast to the case of a regular photonic band edge, in the 
vicinity of loq both Bloch contributions to develop strong singularity. Specifically, 
as the frequency uj approaches cJo f both contributions grow dramatically, while 
remaining nearly equal and opposite in sign at the slab boundary [H] 

^i pr (0) w -*ei; (0) oc \u> - LJo\~ 1/3 , as w w . (25) 

Due to the destructive interference the resulting field 

1-t (0) = * pr (0) + (0) 

at the surface at z = is small enough to satisfy the boundary condition (@J. As the 
distance z from the slab boundary increases, the destructive interference becomes less 
effective - in part because the evanescent contribution decays exponentially 

(z) w V ev (0) exp (-zk") , (26) 

while the amplitude of the propagating contribution remains constant and very 
large. Eventually, the transmitted wave ^>t (z) reaches its large saturation value 
corresponding to its propagating component ^ pr , as seen in Fig. EJa). 
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Note that the imaginary part k" of the evanescent mode wavenumber in l|26[l also 
vanishes in the vicinity of stationary inflection point 

k" oc \lu — cjol 1 ^ 3 , as lu — > ujq, (27) 

reducing the rate of decay of the evanescent contribution H26[) . As a consequence, the 
resulting amplitude (z) of the transmitted wave reaches its large saturation value 
typ r in H25fl only at a certain distance Z from the surface. 

Z oc 1/k" oc \uj - uj \~ 1/3 . (28) 

This characteristic distance diverges as the frequency approaches its critical value uiq. 

If the frequency of the incident wave is exactly equal to the frozen mode frequency 
wo, the transmitted wave (z) does not reduce to the sum (|22|l of propagating 
and evanescent contributions, because at w = ujq, there is no evanescent solutions 
to the Maxwell equations Instead, ^>t (z) corresponds to a non-Bloch Floquet 

eigenmode diverging linearly with z UJ. 

<f T (z) - * T (0) oc z* , atw=w - (29) 
Such a solution is shown in Fig. ffic). 




Figure 9. Destructive interference of the propagating and evanescent components 
of the transmitted wave inside semi-infinite photonic crystal. The frequency is 
close but not equal to that of stationary inflection point, (a) The squared modulus 
of the resulting transmitted field - its amplitude at z = is small enough to 
satisfy the boundary conditions (b) the squared modulus of the propagating 
contribution, which is independent of z\ (c) the squared modulus of the evanescent 
contribution, which decays with the distance z. The amplitude of the incident 
wave is unity. The distance z from the surface is expressed in units of L. 



2.2.3. Degenerate band edge While the situation with a regular photonic band edge 
ifTTJl appears trivial, the case of a degenerate band edge proves to be quite different. 
Just below the degenerate band edge frequency lu^ (inside the transmission band), 
the transmitted field is a superposition (|22J) of one propagating and one evanescent 
components. Above (inside the band gap), the transmitted wave is a combination 
H21fl of two evanescent components. In this respect, a regular and a degenerate band 
edges are similar to each other. A crucial difference, though, is that in the vicinity 
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of a degenerate band edge, both Bloch contributions to the transmitted wave diverge 
as to approaches u>d, both above and below the band edge frequency. This constitutes 
the frozen mode regime. 

Let us start with the transmission band. As the frequency lj approaches L04 — 0, 
both Bloch contributions in ill'l'l) grow sharply, while remaining nearly equal and 
opposite in sign at the surface 

* pr (0) M -* e „(0) OC |w d -^r 1/4 , HSLU^LU d -0. (30) 

This asymptotic formula was obtained in |13) using the perturbation theory for the 
4x4 transfer matrix (|67[) . The destructive interference (|30|l ensures that the boundary 
condition J3J can be satisfied, while both Bloch contributions to vpy (z) diverge. As 
the distance z from the slab boundary increases, the evanescent component ^ ev (z) 
dies out 

# e „ (z) « * e „ (0) exp {-zk") (31) 

while the propagating component typr (z) remains constant and very large. Eventually, 
as the distance z further increases, the transmitted wave ^t (z) reaches its large 
saturation value corresponding to its propagating component ^ pr (z), as illustrated in 
Fig. EH Note that the imaginary part k" of the evanescent mode wavenumber also 
vanishes in the vicinity of degenerate band edge 

k" oc \oj — uJd^^ j as lu — > LUd, (32) 

reducing the rate of decay of the evanescent contribution l|31|l . As a consequence, the 
resulting amplitude "Jt {%) of the transmitted wave reaches its large saturation value 
typ r only at a certain distance Z from the surface 

Z cx 1/k" oc \w -uj d \~ l,i . (33) 

This characteristic distance increases as the frequency approaches its critical value u>d, 
as illustrated in Fig. 0Ja) and (b). 

If the frequency uj of the incident wave is exactly equal to u)d, the transmitted wave 
(z) does not reduce to the sum of two Bloch contributions. Instead, it corresponds 
to a non-Bloch Floquet eigenmode linearly diverging with z 

* T (z) - * T (0) oc z^ d , at ui = Lu d - (34) 

This situation is shown in Fig. life). 

The above behavior appears to be very similar to that of the frozen mode regime 
at a stationary inflection point, shown in Figs. and Yet, there is a crucial 
difference between the frozen mode regime at a stationary inflection point and at a 
degenerate band edge. According to (|19J) . in the immediate proximity of a degenerate 
band edge, the Pointing vector St of the transmitted wave is infinitesimal, in spite 
of the diverging wave amplitude (|30|) . In other words, although the energy density 
Wt oc \^t\ 2 of the frozen mode diverges as u> — > uid — 0, it does not grow fast enough 
to offset the vanishing group velocity 1(17)1. As a consequence, the photonic crystal 
becomes totally reflective at oj = uid- Of course, the total reflectivity persists at 
bj > uid, where there is no propagating modes at all. By contrast, in the case (I2"§1) of a 
stationary inflection point, the respective Pointing vector St is finite and can be even 
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close to that of the incident wave, implying low reflectivity and nearly total conversion 
of the incident wave energy into the frozen mode. 

The character of frozen mode regime is different when we approach the degenerate 
band edge frequency from the band gap. In such a case, the transmitted field "fy (z) 
is a superposition (|?TJl of two evanescent components. As the frequency w approaches 
cod, both evanescent contributions grow sharply, while remaining nearly equal and 
opposite in sign at the photonic crystal boundary 

^evi (0) « -* eu2 (0) oc l^-wf 174 , asw^w d + 0. (35) 

This asymptotic formula also was derived using the perturbation theory for the 
4x4 transfer matrix l|67[) . Again, the destructive interference (I35|l ensures that the 
boundary condition (0J is satisfied, while both evanescent contributions to {z) 
diverge in accordance with (|35J) . As the distance z from the slab boundary increases, 
the destructive interference of these two evanescent components is lifted and the 
resulting field amplitude increases sharply, as seen in Fig. Illf a). But eventually, 
as the distance z further increases, the transmitted wave (z) completely decays, 
because both Bloch contributions to (z) are evanescent. The latter constitutes the 
major difference between the frozen mode regime above and below the DBE frequency 
0Jd- The rate of the amplitude decay, as well as the position of the maximum of the 
transmitted wave amplitude in Figs. 0Jd - f ) and lllf al , are characterized by the 
distance Z in (|33fl . 




20 40 60 20 40 60 20 40 60 

Distance z Distance z Distance z 

Figure 10. Destructive interference of the two Bloch components of the 
transmitted wave inside semi-infinite photonic crystal. The frequency is w = 
ujd — 10 _4 c/L, which is slightly below the degenerate band edge in Fig. I15f b~), 
(a) The squared modulus of the resulting transmitted field - its amplitude at z = 
is small enough to satisfy the boundary conditions Jl}; (b) the squared modulus of 
the propagating contribution, which is independent of z; (c) the squared modulus 
of the evanescent contribution, which decays with the distance z. The amplitude 
of the incident wave is unity. Similar graphs related to the stationary inflection 
point are shown in FiglSl 



2.2-4- Physical reason for the growing wave amplitude If the frequency lu is close, 
but not equal, to that of a stationary point © of the dispersion relation, the wave 
^>t (z) transmitted to the semi-infinite periodic layered medium is a superposition of 
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Figure 11. Destructive interference of the two Bloch components of the 
transmitted wave inside semi-infinite photonic crystal. The frequency is w = 
u>d + 10 -5 c/L, which is just above the degenerate band edge in Fig. I15f b). (a) 
The squared modulus of the resulting transmitted field - its amplitude at z = 
is small enough to satisfy the boundary conditions (b) and (c) the squared 
moduli of the two evanescent contributions; both decay with the distance z. The 
amplitude of the incident wave is unity. Similar graphs related to the frequency 
just below DBE in Fig. ITSTbl are shown in Fig. I1UI 



two forward Bloch modes \&i (z) and \& 2 (z) 

* T (z) = #i (z) + * 2 (z) . (36) 

The two Bloch modes in (|36|l can be a propagating and an evanescent, as in l|22|l. 
or they can be both evanescent, as in (|21|l . In the vicinity of frozen mode regime, 
as the frequency approaches its critical value (u>o or ujd), the two Bloch eigenmodes 
contributing to ^>t (z) become nearly indistinguishable from each other 

*i (z) w a* 2 (z) , asw^ oj s , (37) 

where a is a scalar, and oj s is the frozen mode frequency (uiq or oj^). The asymptotic 
relation (|37|l reflects a basic property of the transfer matrix (|61|l of periodic layered 
structures at frequency of either a stationary inflection point, or a degenerate band 
edge. A rigorous analysis based on the perturbation theory and leading to (|37|l is 
carried out in ^21- For more on this see the Section 4. 

Let us show under what circumstances the property l|37(l can lead to the frozen 
mode regime. The sum (|36|) of two nearly parallel column vectors "Ji and ^2 must 
match the boundary conditions (@J with the incident and reflected waves. If the 
incident wave polarization is general, then the nearly parallel Bloch components ^1 
and ^2 must be very large and nearly equal and opposite 

*i (0) « -* 2 (0) , |*i (0)| « |*2 (0)| » |*/|, (38) 

in order to satisfy the boundary conditions Indeed, since the incident field 

polarization is general, we have no reason to expect that the column vector ^ (0) 
at the surface is nearly parallel to *i (0) and *2 (0). But on the other hand, the 
boundary conditions say that 



* (0) = (0) + * 2 (0) 



(39) 
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Obviously, the only situation where the sum l|39fl of two nearly parallel vectors can be 
not nearly parallel to either of them is the one described in (|38l) . 

There is one exception, though. As we already stated in (|37|l , in the vicinity of the 
frozen mode frequency, the two Bloch components "Ji and \&2 of the transmitted wave 
are nearly parallel to each other. For this reason, if the polarization of the incident 
wave ^1 is such that ^ (0) in 1)39(1 is nearly parallel to one of the Bloch eigenmodes 
\&i (0) or ^2 (0), it is also nearly parallel to the other one. So, all three column vectors 
*i (0), ^2 (0), and ^ (0) are now parallel to each other. In this, and only this case, the 
amplitude of the transmitted wave ^>t (z) will be comparable to that of the incident 
wave. There is no frozen mode regime for the respective vanishingly small range of 
the incident wave polarization. A particular case of the above situation is a regime of 
a single mode excitation, where only one of the two Bloch components ^\ or $2 i n 
l(36J) contributes to the transmitted wave. 

Finally, let us reiterate that in the limiting cases of lo = loq or u> = ujd, the 
transmitted wave (z) corresponds to the non-Bloch Floquet eigenmode l|29|) or i|34[) . 
respectively. Either of them linearly diverges with z. Again, the only exception is when 
the incident wave has the unique polarization, at which the transmitted wave (z) 
is a propagating Bloch eigenmode with zero group velocity and a limited amplitude, 
comparable to that of the incident wave. Incident wave with any other polarization 
will generate the frozen mode inside the periodic medium. 

2. 3. Frozen mode regime in bounded photonic crystals 

The above consideration was based on the assumption that the transmitted wave is 
a superposition of only forward waves, which include propagating modes with u > 
and/or evanescent modes with k" > 0. This assumption, reflected in l|2UI) . (|21|l . and 
<|22[1 . does not apply to bounded photonic crystals, where the periodic medium does 
not occupy the entire half-space z > 0. If the periodic layered array in Fig. ^has a 
finite thickness, the electromagnetic field inside the periodic stack is a superposition 
of all four Bloch eigenmodes with either sign of the group velocity u of propagating 
contributions, and/or cither sign of k" of evanescent contributions. How does a finite 
thickness affect the frozen mode regime? In Fig. \5\ we depicted the frozen mode profile 
in periodic stacks composed of different number N of identical unit cells L in Fig. H3 
The total thickness of the respective photonic slab is equal to NL. In all cases, the 
incident wave frequency u> coincides with that of the degenerate band edge in l|ll(l . 
The presence of the second (right-hand) boundary of the periodic array gives rise to 
the backward wave contribution to {z). Comparison of Fig. [SJa - c) to the semi- 
infinite case in Fig. |5{d) shows that the backward wave contribution to the formation 
of the frozen mode profile becomes significant only at a certain distance from the 
surface of incidence at z = 0. Specifically, the backward wave contribution eliminates 
the frozen mode in the right-hand portion of the finite photonic slab, while having no 
impact in its left-hand portion at < z <C NL. Similar situation occurs at the frozen 
mode regime associated with stationary inflection point (Q. Importantly, the frozen 
mode profile near the surface of incidence is not affected by the finite dimensions of 
the photonic crystal. 

In addition to the modification of the frozen mode profile, the bounded 
photonic crystals of relatively small dimensions can display strong Fabry-Perot cavity 
resonances, also known as transmission band edge resonances. The respective 
resonance frequencies are located strictly inside the transmission band (see, for 
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example, |14j . and references therein) and do not interfere with the frozen mode regime. 
Cavity resonances are distinct from the frozen mode regime and go outside the scope 
of our investigation. 

3. Frozen mode regime at oblique propagation abnormal grazing modes 

A phenomenon similar to the frozen mode regime can also occur at oblique wave 
propagation, where the incident, reflected and transmitted waves are all propagate at 
an angle to the z axis, as shown in Fig. 1121 Consider the situation where the normal 
component u z of the group velocity of the transmitted propagating wave vanishes, 
while the tangential component uj_ remains finite. 

u z = = 0, u± = ^ 0, at uj = u> s = to (ks) . (40) 
ok z dk± v > 

This is exactly what happens in the vicinity of the well-known phenomenon of total 
internal reflection [jjj. Similar effect occurs in any photonic crystal at frequency 
corresponding to the transmission band edge for a particular direction of incidence. 
Remarkably, the total reflection of the incident wave is not the only possible outcome. 
Another alternative is that the transmitted wave forms an abnormal grazing mode 
with dramatically enhanced amplitude and tangential energy flux. The profile of such 
a grazing mode, i.e., the field dependence on the distance z from the surface, appears 
to be very similar to that of the frozen mode at normal incidence shown in Figs. |5] 
and 0] The only difference is that the tangential component of the transmitted wave 
group velocity now remains finite and can even be comparable to the speed of light in 
vacuum. 

A significant advantage of the oblique version of the frozen mode regime is that it 
can occur in much simpler periodic structures, compared to those supporting the frozen 
mode regime at normal incidence. Examples of periodic layered arrays supporting 
only the oblique version of the frozen mode regime are shown in Figs. and |H1 
These structures are too simple to support any kind of frozen mode regime at normal 
incidence - they have only two different layers in a unit cell L, of which only one layer 
is anisotropic. But at oblique incidence, these relatively simple periodic arrays can 
support the frozen mode regime. The presence of at least one anisotropic layer in a 
unit cell L is still required. The physical requirements to periodic structures capable of 
supporting both normal and oblique versions of the frozen mode regime are discussed 
in Section 5. 

3.1. Axial dispersion relation - basic definitions 

Consider a plane monochromatic wave obliquely incident on a periodic semi-infinite 
stack, as shown in Fig. El Due to the boundary conditions J3J, the incident, reflected, 
and transmitted waves should be assigned the same pair of tangential components 
k x , k y of the respective wave vectors 

fc). = fa). - fa)„. fr), = H ^ H ' (41) 

while the axial (normal) components k z are all different. For the incident and reflected 
waves we have simply 



(42) 
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Figure 12. Scattering problem for a plane wave obliquely incident on a semi- 
infinite periodic layered medium. The arrows schematically shows the Pointing 
vectors of the incident, reflected and transmitted waves. The amplitude of the 
incident wave is unity. 
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Figure 13. The case of oblique incidence, where the transmitted wave is a grazing 
mode with tangential energy flux and a frozen mode profile. 
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Let us turn to the transmitted wave. The transmitted wave is usually a 
composition of two Bloch eigenmodes with the same k± = (k x ,k y ) from 1)41(1 . but 
different k z and different polarizations. For given k± and lo, the value of k z is obtained 
by solving the time-harmonic Maxwell equations (|52ll in the periodic medium. The 
so-obtained correspondence between the wavenumber k z and the frequency lo at fixed 
fcj. is referred to as the axial or normal dispersion relation. Real k z correspond 
to propagating (traveling) Bloch modes, while complex k z correspond to evanescent 
modes, decaying with the distance z from the surface. Unlike k x and k y , the Bloch 
wavenumber k z is defined up to a multiple of 2t:/L. 

Similarly to the case of normal propagation, the expression (|40|l defines stationary 
points of the axial dispersion relation. The definition l|4U|l is a generalization of (JSJ) to 
the case of oblique propagation. Different kinds of axial stationary points are defined 
as follows. 

- A regular band edge of axial dispersion relation 

- A stationary inflection point of axial dispersion relation 

doj d 2 oj d 3 uj . . 

Qk z =^ W t = °' dkl^°- (44) 

- A degenerate band edge of axial dispersion relation 

8lo 3 2 lu 3 lo d 4 oj , 

wr^m = ^M =Q >M^ (45) 

The above definitions are analogous to those in 1(11(1 . (JJJ, and (J3J), related to 
the case of normal propagation. We still can refer to the band diagrams in Fig. 
where the quantity k is now understood as the normal component k z of the Bloch 
wavenumber at fixed k±. 



3.2. Grazing mode solutions 

All basic features of axially frozen mode regime at oblique incidence are virtually the 
same as in the case of normal propagation. In particular, all the expressions (12011 
though ((39(1 of the previous section describing the structure and composition of the 
transmitted field (z) remain unchanged. This close similarity holds both for the 
frozen mode regime at a stationary inflection point (|44|l and at a degenerate band edge 
H45fl . In either case, Figs. and 0] give an adequate idea of the frozen mode profile. 
Still, there is one essential qualitative difference. Namely, in the case of axially frozen 
mode we have to remember that the tangential component of the group velocity is not 
zero, even if the normal component 1)40(1 vanishes. This means that the axially frozen 
mode is in fact an abnormal grazing mode with purely tangential energy flux, greatly 
enhanced amplitude, and a very unusual profile, similar to that shown in Figs. [21 and 
0] The steady-state tangential energy flux of the transmitted wave is 

S x =W T (z)u x , (46) 

where the tangential component ui of the group velocity remains large in the vicinity 
of axially frozen mode regime. Therefore, the tangential energy flux S± also grows 
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dramatically, as the frequency approaches the respective critical point l|44|l or 1451) . 
This situation is illustrated in Fig. ED 

Note that the purely tangential energy flux in the transmitted wave in Fig. If 3l does 
not mean that this grazing mode can be classified as a surface wave. Indeed, a surface 
wave is supposed to decay with the distance from the interface in either direction, 
which is not the case here. The possibility of abnormal surface waves associated 
with a degenerate band edge of axial dispersion relation will be addressed in the next 
subsection. 

3.3. Subsurface waves in the vicinity of degenerate band edge of axial dispersion 
relation 

So far in this section we have tacitly assumed that 

k 2 x +k 2 y <u> 2 c 2 . (47) 

The inequality H47JI implies that the z component (|42J) of the wave vector of the 
incident wave is real. This is a natural assumption when considering the problem of 
a plane wave incident on a semi-infinite photonic crystal. 
Consider now the opposite case where 

k 2 x + k 2 y > oj 2 c 2 . (48) 

In this situation, there is no plane propagating waves in vacuum matching the 
boundary conditions (|4*T|) . Still, if the frequency u lies inside a band gap for a given 
there can be a solution for ty? corresponding to a surface wave (see, for example, 
P"5*| and references therein). Generally, such a solution is a superposition l|21|) of two 
evanescent modes. 

Consider now a surface wave at frequency located inside a photonic band gap and 
close to a degenerate band edge ()45|l of the axial dispersion relation for a given k± . 
The Bloch composition of the field inside the periodic medium is 

* T (z) = Vevi (z) + ^ev2 (z) , where u>u d , z > 0. (49) 

As frequency lo approaches uJd, the column vectors >t e „i and ^ e v2 in l|49f) become 
nearly parallel to each other (see (|37|l and comments therein). Together, they can 
form a surface wave with the profile similar to that of the frozen mode shown in Fig. 
Illf a). Although formally, it would still be a surface wave, its profile is highly unusual. 
Namely, the field amplitude inside the periodic medium sharply increases with the 
distance z from the surface, reaches its maximum at a certain distance Z defined in 
(|33|l . and only after that it begins a slow decay. Since the field amplitude reaches its 
maximum only at a distance from the surface, and the respective maximum value can 
exceed the field amplitude at the interface by several orders of magnitude, such a wave 
can be referred to as a subsurface wave. 

4. Floquet modes at stationary points of dispersion relation 

Whether or not a given photonic crystal can support the frozen mode regime is 
determined by its (axial) electromagnetic dispersion relation. Specifically, if the 
dispersion relation develops a stationary inflection point or a degenerate band edge, 
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then one can always expect the frozen mode regime in the vicinity of the respective 
frequency. Restricting ourselves to periodic layered structures, we can link the 
symmetry of the periodic array to the possibility of the existence of the proper 
stationary point of the dispersion relation. 

The first subsection of this section starts with some basic definitions and 
notations of electrodynamics of stratified media involving birefringent layers. We 
briefly describe the formalism of 4 x 4 transfer matrix, generalized to the case of 
oblique wave propagation. Different modifications of this approach have been used 
in electrodynamics of stratified media for at least two decades (see, for example, ^5] 
and references therein). Wherever possible, we use exactly the same notations and 
terminology as in |51ll0|. 

In the second subsection we establish the relation between the symmetry of the 
periodic layered array and the possibility of the existence of a degenerate band edge in 
the respective dispersion relation. The emphasis is on the case of oblique propagation, 
where the symmetry restrictions on the periodic array are much less severe. As a 
consequence, the axial frozen mode regime at oblique incidence can occur in periodic 
structures that are too simple to support the frozen mode regime at normal incidence. 
Examples of the periodic layered structures supporting the (axial) dispersion relation 
with a degenerate band edge are considered in the next section. 

Similar problem for the case of a stationary inflection point was addressed in 
|10l ITT] . Note that the conditions for the existence of a stationary inflection point and 
a degenerate band edge are mutually exclusive. 

4-1- Time-harmonic Maxwell equations in periodic layered media 
Our analysis is based on time-harmonic Maxwell equations 

V X E (x, y, z) = i—B (x, y,z) , V x H (x, y, z) = —i—D (x, y, z) , (50) 
c c 

with linear constitutive relations 

D (x, y,z)=e (z) E (x, y,z), B (x, y,z) = p, (z) H (x, y, z) . (51) 

In layered media, the material tensors e and fi in H51(l depend on a single Cartesian 
coordinate z. Using Ij51|) . the Eqs. (|50J) can be recast as follows 

V x E (x, y, z) = i-pL (z) H(x,y,z), V x H (x, y, z) = -i-e (z) E (x, y, z) . (52) 
c c 

Let us turn to the scattering problem of Fig. 1121 Given the boundary conditions 
H41f> . the field dependence in (|52|) on the transverse coordinates x and y can be 
accounted for by the following substitution 

E (r) = e t{k * x+k y^£ (z) , H (r) = e i( ~ k ' x+k ^n (z) , (53) 

which also allows to separate the tangential field components into a closed system of 
four linear differential equations 



d z ^{z) =i-M{z)^{z) 



(54) 
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"3/ (z) in l|54|) is a vector-column 

*(*) 



£r (z) 

Sy(z) 

£ x {z) 

8y{z) 



(55) 



The normal field components E z (z) and H z (z) can be expressed in terms of ^ (z). 
In the particular case of k \\ z, (|55J) turns into 

The system l|54l) is referred to as the reduced Maxwell equations. It is relevant 
only if the time-harmonic electromagnetic held can be assigned a certain value of 
k± = (k x ,k y ), which is case here due to the boundary conditions (|4ip. The 4x4 
matrix M (z) in (|54|l is referred to as the (reduced) Maxwell operator. Note that 
due to the substitution l|53(l . the Maxwell operator M (z) depends not only on the 
physical parameters of the periodic structure, but also on the tangential components 
k± = (k x ,k y ) of the wave vector. The explicit expression for M (z) for the case 
of oblique propagation in stratified media composed of birefringent layers is rather 
cumbersome. It can be found , for example, in [TJ]|, along with extensive discussion 
of its analytical properties. 

4-2. The transfer matrix formalism 
The Cauchy problem 

— *(z)=i-M(z)*{z),H>(zo) = V (56) 
oz c 

for the reduced Maxwell equation l|54|l has a unique solution 

9(z)=T(z,zo)9(zo). (57) 

The 4x4 matrix T (z, Zq) is referred to as the transfer matrix. It relates the values of 
time-harmonic electromagnetic field (z) at any two points zq and z of the stratified 
medium. The transfer matrix of a stack of layers is defined as 

T s = T(D,0), 

where z = and z = D are the stack boundaries. The transfer matrix of an arbitrary 
stack is a sequential product of the transfer matrices T m of the constitutive layers 

T s = \{T m . (58) 

m 

If the individual layers m are uniform, the corresponding single-layer transfer matrices 
T rn can be explicitly expressed in terms of the respective Maxwell operators M m 

T m = exp (iD m M m ) , (59) 

where D m is the thickness of the m-th layer. The explicit expression for the Maxwell 
operator M m of an arbitrary uniform anisotropic layer can be found, for example, in 
Therefore, Eq. (|58|l together with l|59(l provide an explicit analytical expression 
for the transfer matrix T$ of an arbitrary stack of uniform dielectric layers with or 
without anisotropy, for an arbitrary (oblique or normal) direction of propagation. 
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4-2.1. Transfer matrix in periodic layered media In a periodic layered medium, the 
4x4 matrix M(z) in (|54|) is a periodic functions of z 

M(z + L)=M (z) . 

Bloch solutions *$>k {z) of the reduced Maxwell equation (|54fl with the periodic M(z) 
are defined as 

(z + L) = e lkL ^ k (z) . (60) 

In the case of oblique propagation, the quantity k in ([60(1 denotes the z component of 
the Bloch wave vector. 

Introducing the transfer matrix of a unit cell L 

T L =T(L,0), (61) 

we have from (O, lHJUJ) and (JHH 

T L V k (0) = e tkL ^ k (0) . (62) 

Thus, the four eigenvectors 

tfi(0), i= 1,2,3,4. (63) 

of the transfer matrix Tl of a unit cell are uniquely related to the Bloch solutions 
^>k {z) of the reduced Maxwell equation l|54|) . The respective four eigenvalues 

X i = e* k * L ,i = \,2,S,A (64) 

of Tl are the roots of the characteristic polynomial 

det (T L - XI) . (65) 

Each of the four eigenvectors l|63ll corresponds to either propagating or evanescent 
Bloch wave, depending on whether or not the respective Bloch wavenumber fc, from 
H64f> is real. Further in this section we will see that the relation 1|63[) does not apply at 
stationary inflection points of the (axial) dispersion relation. This important exception 
is directly related to the very nature of the frozen mode regime. 

The explicit expressions for the 4x4 transfer matrix 1611) , along with the detailed 
description of its analytical properties can be found, for example, in [101 113| . 

4-2.2. Transfer matrix at stationary points of dispersion relation Although at any 
given frequency to, the reduced Maxwell equation (|54[l has exactly four linearly 
independent solutions, it does not imply that all four of them are Bloch waves as 
defined in (|60[l . Specifically, at frequencies of stationary points l[43[l , 144|) , or 145|) where 
the axial component of the group velocity of some of the propagating modes vanishes, 
some of the four solutions can be algebraically diverging with z and, therefore, cannot 
be classified as Bloch waves. For example, at frequency lu^ of a degenerate photonic 
band edge, the four solutions of Eq. I|54[l include a propagating Bloch mode with zero 
group velocity and three Floquet eigenmodes diverging as z, z 2 , and z 3 , respectively 
(see the details in ^^E])- Some of these eigenmodes are essential for understanding 
the frozen mode regime. 

Consider such non-Bloch solutions in terms of the transfer matrix Tl of a unit 
cell. Although the matrix l[61[l is invertible, it is neither Hermitian, nor unitary and, 
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therefore, may not be diagonalizable. Specifically, if the frequency approaches one 
of the stationary points Q43[l. I|44|) . or l|45|) . some of the four eigenvectors (0) m 
16211 become nearly parallel to each other. Eventually, as u> reaches the stationary 
point value, the number of linearly independent eigenvectors (0) becomes smaller 
than four, and the relation l|63[l does not apply at this particular point. The number 
of linearly independent eigenvectors of the transfer matrix is directly linked to its 
canonical Jordan form. 

At a general frequency different from that of any stationary point of the (axial) 
dispersion relation, the transfer matrix Tl is diagonalizable, and its canonical Jordan 
form is trivial 



a ik2 



Tl(lo) 



o 
o 







(66) 



This matrix has four linearly independent eigenvectors (Mi.'St corresponding to four 
Bloch eigenmodes, each of which is either propagating, or evanescent. The four 
respective values of the wavenumber k are determined by Ij64|l and f55} . 

At frequency u>d of a degenerate band edge, the canonical Jordan form of the 
transfer matrix Tl becomes 



e lk 







1 

e lkd 



(67) 
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e ikd 
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gjfcd 
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associated with one propagating Bloch mode with zero group velocity, and three non- 
Bloch eigenmodes diverging as z, z 2 , and z 3 , respectively . If the frequency uj deviates 
from that of the stationary point, the transfer matrix Tl becomes diagonalizable 
with the canonical Jordan form (|66l) . The perturbation theory relating the non-Bloch 
eigenmodes at the frequency of a degenerate band edge to the Bloch eigenmodes in 
the vicinity of this point is presented in |13| . 

The other possibilities include a regular band edge l|43[l and a stationary inflection 
point Q44J) • Those cases are discussed in ^3 EH • 



5. Periodic layered structures with degenerate band edge of axial 
dispersion relation 

Not any periodic stack can have electromagnetic dispersion relation with a degenerate 
band edge (|43t . One fundamental restrictions stems from the fact that at the frequency 
uid of a degenerate band edge, the transfer matrix Tl must have the canonical form 
(|67|l . Such a matrix cannot be reduced to a block-diagonal form, let alone diagonalized. 
Therefore, a necessary condition for the existence of a degenerate band edge is that 
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Wave number k Wave number k Wave number k 

Figure 14. Dispersion relation u} (k) of the periodic stack in Fig. |H] for three 
different values of the misalignment angle cj>. In the cases <f> = (no misalignment) 
and rf> = 7r/2, none of the spectral branches can develop a degenerate band edge 
(DBE). While in the case (c) of = 7r/4, one of the spectral branches develops a 
DBE. 



the symmetry of the transfer matrix Tl does not impose its reducibility to a block- 
diagonal form. The above condition does not imply that the transfer matrix Tl must 
not be reducible to a block-diagonal form at any frequency w. Indeed, at a general 
frequency w, the matrix Tl is always reducible and even diagonalizable. The strength 
of the symmetry imposed reducibility, though, is that it leaves no room for exceptions, 
such as the frequency aid of degenerate band edge, where the transfer matrix Tl must 
not be reducible to a block-diagonal form. Therefore, in the case of symmetry imposed 
reducibility of the transfer matrix, the very existence of degenerate band edge H45|) is 
ruled out. 

At this point we would like to emphasize the important difference between the 
cases of normal and oblique propagation. In the case of normal propagation, the 
symmetry of the Maxwell operator M (z) in l)54|l and the transfer matrix Tl simply 
reflects the symmetry of the periodic layered array. By contrast, in the case of oblique 
propagation, the substitution (|53|l lowers the symmetry of the matrices M (z) and 
Tl and makes it dependent not only on the periodic structure itself, but also on the 
orientation of fcj_ = (k x ,k y ) in the x — y plane. Lower symmetry of the "oblique" 
transfer matrix Tl may remove its symmetry-imposed reducibility to a block-diagonal 
form, even if at normal incident such a reducibility was imposed by the symmetry 
group of the periodic structure. In other words, even if a certain periodic layered 
structure cannot support a degenerate band edge at normal propagation, it may 
develop such a stationary point at oblique propagation. In this respect, the situation 
with degenerate band edge is reminiscent of that of a stationary inflection point, where 
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Figure 15. k — u> diagram of the periodic stack in Fig. |S]for four different values 
of the B - layer thickness, (a) B/L = 0.711 44, in this case the upper edge of the 
frequency gap develops a DBE. (b) B/L = 0.374 43, in this case the lower edge 
of the frequency gap develops a DBE. (c) B/L = 0.1. (d) B/L = 0, in this case 
the intersecting dispersion curves correspond to the Bloch waves with different 
symmetries - the respective modes are decoupled. 

the cases of normal and oblique propagation are also essentially different from each 
other [SJ^j. Observe, though, that in periodic layered structures, the possibilities of 
a stationary inflection point and a degenerate band edge are mutually exclusive. 

Further in this section we consider specific examples of periodic layered arrays 
supporting a degenerate band edge and the related frozen mode regime. We start with 
the particular case of normal propagation, requiring more complex periodic structures. 

5.1. Degenerate band edge at normal propagation 

The simplest periodic layered structure capable of supporting a degenerate band edge 
at normal propagation is shown in Fig. |SJ A unit cell L contains one isotropic B 
layer and two misaligned anisotropic layers A± and A<i with inplane anisotropy. The 



Frozen light in photonic crystals with degenerate band edge 



30 




2 3 4 

Wave number k 



Figure 16. Axial dispersion relation ui (k z ) of the two-layered periodic stack 
in Fig. |3 The tangential components k x ,k y of the wave vector are fixed at 
k x = k y = 1.9403. The fifth spectral branch develops a degenerate band edge 
1451 at k z = w/L. 



isotropic layers have the thickness B and the dielectric permittivity 



£b 



e B 
s B 
e B 



(68) 



The dielectric permittivity tensors ea in each anisotropic A layer has the form 



^A + Scos2cp 6sin2ip 
S sin 2tp EA — 5cos2(p 

£ 3 



(69) 



where the parameter 5 characterizes the magnitude of inplane anisotropy and the 
angle ip determines the orientation of the anisotropy axes in the x — y plane. All the 
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A layers have the same thickness A and the same magnitude 8 of inplane anisotropy. 
The only difference between the adjacent anisotropic layers Ai and A 2 in Fig. [fj] is 
their orientation if. 

An important characteristic of the periodic structure in Fig. EJis the misalignment 
angle 

4> = <Pi - f2 (70) 
between the layers A\ and A 2 . This angle determines the symmetry of the periodic 
array and, eventually, the kind of k — to diagram it can display. Fig. 1141 illustrates the 
relation between the misalignment angle cf> and the symmetry of the respective k — uj 
diagram. Generally, there are three possibilities, reflected in Table. 1. 

Table 1. 



Value of 4> 


Symmetry class 


Spectral properties 


Example 


cj) = 


mmm = D 2 h 


x and y polarizations are separated 


Fig.maT 


cj) = tt/2 


Amm = Did 


polarization degeneracy 


Fig. 114Tb) 


0^0,^/2 


222 = D 2 


no polarization degeneracy/separation 


Fig. HUc) 



In the case <p = 0, all anisotropic layers have aligned in-plane anisotropy. The 
term " aligned" means that one can choose the directions of the in-plane Cartesian axes 
x and y so that the permittivity tensors in all layers are diagonalized simultaneously. 
In this setting, the Maxwell equations for the electromagnetic waves with x- and the y - 
polarizations propagating along the z axis are uncoupled, implying that the respective 
transfer matrix can be reduced to the block-diagonal form 



Til 


Tl 2 










T22 














T33 


T34 








T43 


T44 



(71) 



The two blocks in l|71|l correspond to the x and y polarization of light. The forth 
degree characteristic polynomial l|65|l of the block-diagonal matrix (|71|l factorizes into 
the product 

F 4 (X) = F x (X)F y (X), (72) 

where F X (X) and F y (X) are independent second degree polynomials related to 
electromagnetic waves with the x- and the y - polarizations, respectively, propagating 
along the z direction. The k — to diagram for this case is shown in Fig. I14f a). where 
each spectral curve relates to a specific linear polarization of light. In this case, the 
symmetry imposed reducibility of the matrix Tj, to a block-diagonal form (|71|l rules 
out the existence of a degenerate band edge. 

In the case <j> = tt/2, the anisotropy axes in the adjacent layers A\ and A 2 are 
perpendicular to each other. The point symmetry group of the periodic array is now 
D 2 d i which is a tetragonal symmetry class. The tetragonal symmetry results in 
polarization degeneracy, implying that the respective transfer matrix can be reduced 
to the following block-diagonal form 



T L = 



Tn 


T12 








T21 


T22 














Tn 


T12 








T21 


T22 



(73) 
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The two identical blocks in l|73|l correspond to either polarization of light. The forth 
degree characteristic polynomial i|65|l of the block-diagonal matrix (|73|l factorizes into 
the product 

F 4 (X) = F 2 {X)F 2 (X), (74) 

where F 2 (X) is a second degree polynomial related to electromagnetic waves with 
either polarization propagating along the z direction. The k — lu diagram for this case 
is shown in Fig. IMf bh where each spectral branch is doubly degenerate with respect 
to polarization. In this case, the symmetry imposed reducibility of the matrix to 
a block-diagonal form (|73*)l also rules out the existence of a degenerate band edge. 

Finally, in the case 4> ^ 0,7r/2, the periodic stack in Fig. 0has a chiral point 
symmetry described as D 2 . There is no symmetry prohibition of a degenerate band 
edge in this case, because the Bloch modes with different polarizations now have 
the same symmetry and, therefore, are coupled. In this case, one can adjust the 
misalignment angle <f> and/or the relative layer thickness 

b = B/L = B/{2A + B), (75) 

so that a given spectral branch develops a degenerate band edge. The respective value 
of the wavenumber is either k = 0, or k = tt/L. 

In the numerical example in Fig. 1151 we show four k — ui diagrams of the periodic 
structure in Fig. corresponding to four different values of the ratio b in 1(751 . In all 
cases, the misalignment angle tfi is equal to 7r/4. The k — uj diagrams in Figs. I15f a) 
and !15f b) show a degenerate band edge in the respective spectral branches. The k — oj 
diagrams in Figs. IT4T c) and lTST bl are identical. 

If the isotropic B layers are completely removed from the periodic structure 
in Fig. El the point symmetry group of the periodic array rises from D 2 to D 2 h, 
acquiring a glide mirror plane my . The two different linear polarizations now become 
uncoupled regardless of the misalignment angle (f>, while the transfer matrix Tl of the 
stack displays a symmetry imposed reducibility to a block-diagonal form l|71|) . The 
respective k — lu diagram is shown in Fig. I15f d^ . 



5.2. Degenerate band edge at oblique propagation 



Consider now a periodic structure with just two layers A and B in a unit cell, as shown 
in Fig. [5] The dielectric material of the A layer has an inplane anisotropy l|69|) while 
the B layer is isotropic. For specificity, we can set 



£a = 



e A + S 
e A - 5 
e 3 



e B 
e B 
e B 



(76) 



Note that the structure in Fig. H3 reduces to that in Fig. [HI in the particular case of 
= 0. Indeed, if the misalignment angle between A\ and A 2 in Fig. El is zero, these 
two anisotropic layers together make a single A layer with double thickness. Therefore, 
at normal propagation, the k — to diagram of the periodic array in Fig. |S]is similar to 
that of the periodic structure in Fig. 0with = 0. The latter is shown in Fig. I14f a^l . 
There is no possibility of a degenerate band edge in this case. 

The situation remains unchanged if the direction of propagation deviates from 
the z axis, but is confined to either the x — z, or the y — z plane. In either case, the 
respective plane is the mirror plane of the transfer matrix, ensuring that the x and 
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y polarizations remain uncoupled. Uncoupled polarizations imply that the transfer 
matrix Tj, is reducible to the block-diagonal form l|71|) at all frequencies. Again, there 
is no possibility of a degenerate band edge in this case either. 

The situation changes only in the case of oblique propagation with k x ,k y ^ 0. 
The Maxwell equations l|54|) for different light polarizations are not decoupled any 
more, and the respective transfer matrix cannot be automatically reduced to a block- 
diagonal form at all frequencies. As a consequence, at certain direction of propagation, 
some spectral branches can develop a degenerate band edge, as shown in the example 
in Fig. 1161 In fact, for any given frequency u>d within a certain frequency range, one 
can find a specific direction k± = (k x ,k y ) for which the degenerate band edge (1451) 
occurs at the chosen frequency u>d- In the example in Fig. 1161 we simply set k x = k y 
and A = B. In such a case, the degenerate band edge frequency and the respective 
value of k x = k y are predetermined by the physical parameters of the periodic array. 



5. 3. Values of physical parameters used in numerical simulations 

In all numerical simulations related to nonmagnetic layered structures in Figs. [S]and 
[H]we use the following values of material parameters in (|68p. (|69|l and l|76|l 



e A = 11-05, S = 7.44, e 3 = 18.49, e B = 1. 



(77) 



The relative thickness of the A and B layers in Fig. as well as the value of the 
misalignment angle (|70|l can be different in different cases. 

Frozen mode profiles presented in Figs. 0] 03 and 1111 are computed for the 
same periodic stack in Fig. with the misalignment angle (j> — tt/A and the ratio 
B/L = 0.37443. The respective k — ui diagram is shown in Figs. I14f c) and ll5f b N |. In 
all cases, the incident wave has unity amplitude and linear polarization with E || y. 
Change in polarization results in the change of the frozen mode amplitude, but it only 
slightly affects the Bloch composition of the frozen mode and its dependence on the 
distance z from the surface of incidence. 

In a single case related to a nonreciprocal periodic layered structure with a 
stationary inflection point (Figs. and EJ) we use the following numerical values 
of the electric permittivity and magnetic permeability tensors of the anisotropic A 
layers and magnetic B layers 
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(78) 



(79) 



The misalignment angle <fi in this case is set to be 7r/4. The respective value of the 
stationary inflection point frequency coq at normal propagation is 0.7515 x c/L. 

In all plots of field distribution inside periodic media at z > we, in fact, plotted 
the following physical quantity 



E{z)-E* (z) + H{z)-H* (z) 



(80) 
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which is the squared field amplitude averaged over a local unit cell. The actual function 
\*f? (z)\ 2 , as well as the electromagnetic energy density distribution W (z), are strongly 
oscillating functions of the coordinate z, with the period of oscillations coinciding 
with the unit cell length L. Given the relation W oc \ty (z)\ 2 , the quantity (|80|l can 
also be qualitatively interpreted as the smoothed energy density distribution, with the 
correction coefficient of the order of unity. 

In all plots, the distance z, the wave number k, and the frequency ui are expressed 
in units of L, L _1 , and cL^ 1 , respectively. 

6. Conclusion 

In this paper we outlined several different manifestations of the frozen mode regime 
in photonic crystals. Although all our numerical examples relate to periodic layered 
structures, in fact, the frozen mode regime is a universal wave phenomenon. Indeed, 
we can talk about different kinds of wave excitations in low-loss periodic media. But as 
soon as the respective Bloch dispersion relation displays a singularity like a stationary 
inflection point l|44Jl or a degenerate band edge 1)45(1 , we have every reason to expect a 
very similar behavior involving the frozen mode regime. In other words, the possibility 
of the frozen mode regime is determined by some fundamental spectral properties of 
the periodic structure, rather than by the physical nature of the linear waves. If 
a periodic array is relatively simple - for instance, a stratified medium with one 
dimensional periodicity - its frequency spectrum may prove to be too simple to support 
the proper spectral singularity. The more complex the periodic structure is, the more 
likely it will be capable of supporting such a phenomenon. For instance, in the case of 
layered arrays we need birefringent layers and, at normal propagation, at least three 
layers in a unit cell. 

Another important question is how robust the frozen mode regime is. For 
instance, what happens if we introduce a small absorption or structural imperfections. 
Of course, these factors suppress the frozen mode amplitude. But in this respect, 
the frozen mode regime is no different from any other coherent or resonance effects it 
periodic structures. This problem can be addressed at any particular frequency range 
by appropriate choice of the constitutive materials. 

Another fundamental restriction relates to the size of the periodic structure. In 
this paper we assumed that the periodic array occupies the entire half-space z > 0. 
A good insight on what happens to the frozen mode in a finite periodic array is 
given by Fig. |SJ These graphs demonstrate that the frozen mode regime in a finite 
periodic array can be as robust as that in an hypothetical semi-infinite structure. The 
optimal number of layers depends on such factors as the absorption characteristics 
of the constitutive materials, the geometrical imperfections of the periodic array, the 
desired degree of field enhancement in the frozen mode, etc. On the other hand, in 
finite (bounded) photonic crystals, some new resonance phenomena can arise, such as 
transmission band edge resonances [21 El E]- These effects, though, are qualitatively 
different from the frozen mode regime and occur at distinctly different frequencies. 
The transmission band edge resonance in the vicinity of a degenerate band edge was 
studied in [T^ . 
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